Confinement of two-dimensional excitons in a non-homogeneous magnetic field 
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The effective Hamiltonian describing the motion of an exciton in an external non-homogeneous 
magnetic field is derived. The magnetic field plays the role of an effective potential for the exciton 
motion, results into an increment of the exciton mass and modifies the exciton kinetic energy oper- 
ator. In contrast to the homogeneous field case, the exciton in a non-homogeneous magnetic field 
can also be trapped in the low field region and the field gradient increases the exciton confinement. 
The trapping energy and wave function of the exciton in a GaAs two-dimensional electron gas for 
specific circular magnetic field configurations are calculated. The results show than excitons can be 
_0 , trapped by non-homogeneous magnetic fields, and that the trapping energy is strongly correlated 

Q-' ■ with the shape and strength of the non-homogeneous magnetic field profile. 
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_-* ' Two dimensional confinement of excitons or atoms by magnetic fields is an important phenomenon in physics, and 

f~^ is expected to bring a significant improvement to lasers and other functional xlevicesu. Numerous works have been 

^-H ' carried out on the study of the exciton properties in applied magnetic fieldsa Lj. However, for the most part, all the 

0^ | attention has been focused on the influence of homogenous fields on the exciton inner properties. Several theoretical 
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experiments show the possibility of trapping and guiding of atoms by means of non-homogeneous 
Recently, numerous papersB have appeared on the properties of charged particles, like e. g. 
electrons, in different non-homogeneous magnetic field profiles,. For example, quantum mechanical bound states were 
found which had no classical analogue. Christianen et q£3a performed photoluminescense (PL) measurements on 
excitons in the presence of non-homogeneous magnetic fields. Strips of ferromagnetic material (see, e.g., RefJla) on 
top of a quantum well were used to create strong magnetic field gradients. They concluded that excitons are forced 
to regions of low field gradient and that magnetic traps for excitons are feasible. To the best of our knowledge, there 
is no theoretical work in the literature on studies of excitons in non- homogeneous magnetic fields. 

The purpose of the present paper is to develop a mathematical formalism for calculating the quantum mechanical 
properties of excitons in a non-homogeneous magnetic field, and to illustrate the trapping possibilities of the excitons 
in some magnetic field profiles. These ones are created by the deposition of magnetic disks and also superconducting 
disks on top of a two dimensional electron gas (2DEG) with a homogeneous field applied pecaendicular to the 2DEG. 
It results in a non-homogeneous magnetic dipole type of profile and in a magnetic antidofilj profile, respectively, in 
the 2DEG. 

The layout of the paper is as follows. In Sec. II the effective Hamiltonian describing the exciton motion in the 
non-homogeneous magnetic field is derived. In Sec. Ill we consider circular symmetric magnetic field profiles and 
reduce the exciton Hamiltonian to a one-dimensional (ID) Schrodinger equation subjected to a spatially dependent 
effective potential and effective mass. In Sec. IV we calculate explicitly the magnetic field profiles used in our work 
and study how strongly they are able to trap excitons. Our numerical results for the exciton trapping energies are 
discussed in Sec. V, and the conclusions are given in Sec. VI. 

II. THE EFFECTIVE HAMILTONIAN 

We consider a 2D system of two particles with opposite charge interacting via the Coulomb interaction and moving 
in a non- homogeneous magnetic field characterized by the vector potential A(r). The Hamiltonian of that system is 
given by: 



where e is the dielectric constant of the material the exciton is moving in, m e (rrih) is the electron (hole) effective mass, 
and r e , rh are the electron and hole coordinates, respectively, in the xy-plane. To simplify the above Hamiltonian, we 
use the center of mass R = (m e r e + rrih^hj/M and relative motion coordinates r = r e — r/j, where M = m e + rah is the 
total mass of the exciton. Next, in order to obtain the correct asymptotic .behavior, we apply a wave function phase 
transformation analogous to the one used by Gor'kov and Dzyaloshinskycl in the case of a homogeneous magnetic 
field: 

*(R, r) -► exp {-i(e/hc)r ■ A (R)} *(R, r), (2) 

which leads to the following transformed Hamiltonian: 
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n r . nth . e /_ m e \ emi 



2m h 

„2 



{-^- + 'V,- - iA (a- ^) - >»), + ^a(r,}- 



(3) 

In real experimental situationscIO the size of the exciton is smaller than the length scale over which the magnetic 
field varies, which allows us to use the adiabatic approacha. We assume that the exciton relative motion is fast 
as compared with its center-of-mass motion, and that the characteristic dimension in the relative coordinate (i.e. 
the exciton radius) is much smaller than that of the center-of-mass motion (i.e. the characteristic length of the 
magnetic field inhomogeneity) . This allows us to expand the vector potential A(R± m h i e \v/M) into r power series 
restricting the consideration with up to second order terms. Note that the exciton radius in GaAs. e.g., is typically 
a* B = eh //j,e 2 ~ 120 A, while the magnetic field inhomogeneity varies typically on a micron scaleOE3o. Within the 
adiabatic approximation, we obtain the following expression for the exciton Hamiltonian: 

^=g(-,V r+ ^{(r.V,)A(R) } ) 2 + |l( i V fl + ^[rxB(R)]) 2 -^, (4) 

where /i = m e rrih/M is the exciton reduced mass and £ = (mh — m e )/M. The first term in the Hamiltonian, (Eq. H), 
describes the kinetic energy of the exciton relative motion, in which the expression £(r • Vjj)A(R) can be interpreted 
as an effective vector potential describing the local magnetic field £B(R). We can choose the gauge such that 
(r-Vii)A(R) = B(R) x r/2 + V r {(r • V#) {r • A(R)}} /2, and simplify the second term of this expression by applying 
the following transformation: 

*(R, r) -» exp {-if (e/&c)n(r)} *(R, r), (5) 

with O(r) = (r ■ V^){r ■ A(R)}/2. The above transformation does not change the center-of-mass motion, and leads 
only to the appearance of r 3 -order terms in the relative motion, which are neglected in the present adiabatic approach. 
In doing so, we obtain the following transformed Hamiltonian: 

where terms of first and second power in the magnetic field strength are denoted by 

Wl = 2^ B(R) ' L + iS~c {B(R) X Vr ~ Vr X B(R)} ■ r ' 

W 2 = -^B(R)V. (7) 

Here the symbol L = {rx (— iKV r )} stands for the exciton relative angular momentum operator. Notice that in 
case of a homogeneous magnetic field the above Hamiltonian reduces to the one which can be found in the literature 
(see Ref. 5 and references therein). The exciton Hamiltonian in a homogeneous magnetic field consists of a part with 
no magnetic field dependence, a part linear (the angular momentum term) in the magnetic field, and a quadratic part 



(the diamagnetic shift term) in the magnetic field. In the present Hamiltonian which describes the exciton motion in 
a non-homogeneous magnetic field we also have these terms, but now they are position dependent. We also have an 
additional term with a linear dependence in the magnetic field (the last term in the W\ Hamiltonian, Eq. (M)) which 
results from the gradient in the magnetic field, and which modifies the exciton mass. 

Now within the spirit of the above adiabatic approach, the total exciton wave function is represented as the product 
of the wave functions describing its relative and center-of-mass motions, i. e., ^(R, r) = $(r)i/>(R). Next, it is 
supposed that the relative motion wave function $(r) obeys the Schrodinger equation with the Hamiltonian H rel 
which includes all terms of total Hamiltonian, Eq. (g), except the last one describing the kinetic energy of the 
center-of-mass motion: 

{H rel -E(R : X7 R )}$(r) = 0. (8) 

Note that in contrary to the standard adiabatic approachE-3 the relative motion Hamiltonian (see Eq. (Q)) depends 
not only on the center-of-motion coordinate R but also on the gradient V_r. It does not complicates, however, the 
effective Hamiltonian derivation, but one should take into account that now the eigenvalue of the relative motion 
equation, which usually plays a role of the effective potential for the center-of-mass equation, is an operator and gives 
the correction to the kinetic center-of-mass operator as well. So, the effective center-of-mass Hamiltonian has to be 
presented as follows 
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^ 2M Vfl + jB(R ' Vfl) - (9) 

We solve the relative motion, Eq. (8), by means of a perturbation technique in magnetic field strength powers 
restricting our consideration by the second order terms. This assumption is valid in weak field regime (huj* < 2R*), 

where R* = fj,e 4 /2e 2 h is the effective Rydberg, and ui* — eB/fic is the cyclotron-resanancairequency. This is the 
magnetic field regime relevant for experiments with non-homogeneous magnetic fieldglSLSO which contrasts with 
experiments in homogeneous magnetic fields where usually the high field regime Tito* ^§> 2R* is reached. Notice that 
hcu* = 2R* corresponds to B rs 5 T for GaAs. 

Now taking into account the cylindrical symmetry of the zero order Hamiltonian 

H = ~V 2 r --, (10) 
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we find the zero order wave function 



corresponding to the zero order eigenvalue 

E (n,m) = -X 2 n R* y , A„ = (n - 1/2)" 1 , (12) 

where the symbols n and m stands for the radial and angular relative motion quantum numbers, a* B is the effective 
Bohr radius, A n ^ m is_the normalization constant determined by ($Q ,m (r)|$Q ,m (r)) = 1, and Ln'(x) is the generalized 
Laguerre polynomialEfl. Solving in addition the first order equation 

{H - E Q (n,m)}^ m (r) ^ {E^m) -W 1 }^' m (r), (13) 

with the orthogonality condition ($™' m (r)|$Q' m (r)) = 0, we obtain the first and the second eigenvalue corrections 

E l {n,m) = (*?' TO (r)|WW m (r)> = ^B z (R)m, (14) 

E 2 (n,m) = ($^ m (r)|W«' m (r)> + ($o' m (r)|[^i - £i(n,m)]|$r>)) 

= ^^J^ (R)2 + < 2lP& V « W R ) 2V «) ■ ( 15 ) 

Here the symbols a 7 ^ and /3^ are numerical constants which follow from the averages in Eq. (|l^). For a magnetic 
field perpendicular to the xy-pl&ne we obtained the following values: aj = 21/128, 0q — 3/8, «q = 365/64, [3 2 = 35/4, 
and a 2 ±1 = 145/128, f3 2 ±1 = 11/4. 



Inserting the obtained eigenvalue corrections into Eq. (g) we obtain the final expression for the effective Hamiltonian 

R CM = _il V f 1 _ o»r f^L Bz {nf 
2M H \ mr Mc 2 R; y ' 

- % r K + ^m r B z (R) + ^ r ^-B z (K) 2 , (16) 

which describes the center-of-mass motion of the exciton with relative motion quantum numbers (n r ,m r ) in the 
non-homogeneous magnetic fields. 

III. CYLINDRICAL SYMMETRY 

From now on we limit ourselves to cylindrical symmetric magnetic field profiles B Z (R) = B Z (R). Then, the exciton 
center-of-mass wave function ip(R) can be written as: 

tp(R) = —== exp {im R (f>} ip(R), (17) 

V 2n 

where m# is the quantum number corresponding to the exciton center-of-mass angular momentum. Hamiltonian Eq. 
( [lq ) then leads to the following radial Schrodinger equation: 
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where 
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V en (R)~E}^(R) = 0, (18) 



1 - oC r ciB z (R) 2 
is the effective mass for the center-of-mass motion of the exciton, and 



Meff (R) = 1 .J* „,-»„ , ( 19 ) 
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c i = T^ ; c 2 = — -f- ; <* = —£. (21) 



V e "( R) = (3% r c 2 B z {Rf + m r c 3 B z (R) + E (n r , m r ) + 2Meff , R) -^, (20) 

is the effective potential for the exciton with 

e 2 a* B 2 
R* y Mc 2 "* 8fic 2 ' ^ 2 A ic 

For GaAs we have typically a* B ~ 120 1, R* ~ 5 meV, c x ~ 6.25 x 10" 3 T" 2 , c 2 ~ 0.44 meV T~ 2 and c 3 ~ 0.69 
meV T _1 . Please note that the 2D exciton ground state energy -Eo(l, 0) is — AR* ~ —20 meV and the 2D exciton size 
is one half the Bohr radius, a* B /2 = 60 A. From the above equations, we notice that the non-homogeneous magnetic 
field modifies the exciton center-of-mass motion in the following ways. First the exciton feels a B Z (R) 2 effective 
potential. It implies that excitons will be collected in a region where the magnetic field strength is minimum. Next, 
the exciton mass is enhanced, it increases with B Z (R) 2 , which favours the localization of excitons. Third, the kinetic 
term gives a contribution proportional to the gradient in the _B-field, i.e., a dB z (R)/dR term. Finally, we see that 
the non-homogeneous magnetic field also interacts with the exciton angular momentum, and that this interaction is 
controled by the difference in mass £. The reason is that the exciton is a neutral particle, and if the electron and hole 
would have the same mass, the angular momentum term would not contribute because, electron and hole give the 
same contribution but with opposite sign. 

IV. MAGNETIC FIELD PROFILES 

We solve Eq. (|1S|) numerically using the confinement potential generated by two different non-homogeneous mag- 
netic field profiles, which experimentally are created by the deposition of (a) a magnetized disk and (b) a supercon- 
ducting disk on top of a 2DEG, with a homogeneous magnetic field (B a ) applied perpendicular to the 2DEG. Such 
an experimental configuration results in a non-homogeneous magnetic field profile in the 2DEG. Tha-first profile (a) 
results into a magnetic dipole type of profile, and the second (b) is also called the magnetic antidoiO. A sketch of 
the two experimental systems, the magnetized disk and the superconducting disk, are shown in the inset of Fig. 1(b) 
and Fig. 3(a), respectively. 



A. Magnetized Disk 

To calculate the magnetic field created by a magnetized disk, we assume that the disk is very thin and is completely 
magnetized in the z-direction. Therefore, we can write the magnetization in the following way: 



M(R) = hMS(Z - d)6(a - R)e z , 



(22) 



where h is the disk thickness, a is the disk radius, d the distance of the magnetic disk to the 2DEG, M. the magnetisa- 
tion, 9(a—R) is the Heaviside step function, and Z and R = y/X 2 + Y 2 are cylindrical coordinates. The ccpsesponding 
vector potential can be calculated from the differential form of Amperes law V^A(R) = — 4nV r x M(R)li3. In cylin- 
drical coordinates, the vector potential has only an angular component: 
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where Bq = hM., and K(x) (E(x)) is the elliptic integral of first (second) type. The magnetic field can be evaluated 
straightforwardly from B(R) =Vx A(R), which results into: 



B Z (R) 



(a + R) 



(a + R) 2 + d 2 



A V (R) 



BE- 



R 2 



-|^)-IaV) 



R 2 VaR 

,2\ Q 



l -\)w- K{f?) 



(25) 



The magnetic field profile of the magnetized disk (the magnetic dipole profile) in the 2DEG for d = 0.2 /im, a = 2 //m, 
and Bq — 0.1 T is shown in the inset of Fig. 1(a). Due to the quadratic dependence of the effective potential on the 
magnetic field (see Eq. (22)), we can maximize the confinement potential by applying a homogeneous field B a to the 
magnetic dipole profile, B t ° tal {R) = B a + B Z (R). We use the experimental results of Dubonos et aL3 in order to take 
values for the applied, field (B a ) smaller than the coercivity field, which has a strong dependence with the radius of 
the magnetized diskli3. 

The magnetic dipole profile in the 2DEG for different values of the relation d/a are shown in Fig. 1(a) (Fig. 2(a)) 
as a function of the radial coordinate, in the presence of a homogeneous magnetic field B a = 0.35 T (B a = — 0.25T), 
where we took B§ = 0.1 T. We also show in Fig. 1 (B a = 0.35 T) and in Fig. 2 (B a = -0.25 T) the effective potential, 
Eq. (P0[), for the exciton center-of-mass motion for the Is (Fig. 1(a) and Fig. 2(a)), and for the 2p~ (Fig. 1(b) and 
Fig. 2(b)) exciton relative motion quantum state, and the effective mass (Fig. 1(b) and Fig. 2(b)), Eq. (Jl^), for the 
2p _ exciton relative motion quantum state. 

Notice that the effective potential, Figs. 1-2, is negative. This is due to the energy shift resulting from the zero field 
exciton relative motion Eo(n r , m r ) in the equation of the effective potential (see Eq. (|2C|)). This energy increases with 
increasing principal relative motion quantum number n r (see Eq. (12)). Also notice that the effective potential of the 
excited levels of the exciton relative motion (e.g., see Fig. 1-2 (b) for the 2p~ state) has a confinement region larger 
than the Is relative quantum level. The reason for this is as follows: the coefficient of the dia mag netic term in the 
effective potential equation (the term with quadratic dependence in the magnetic field in Eq. (EQ)), /3^ r , is related 



~(r)|rW») (see 



to the averaging of r 2 in the wave function of the exciton relative motion, i.e., /3^, r oc ($£ 
Eq. (Jlq)) which increases for increasing n r . Then, the confinement for the n r exciton relative motion states should 
be stronger than the one for the n r — 1 relative state. Further, the exciton angular momentum term in the exciton 
Hamiltonian is, in practice, bigger than the diamagnetic term. This is because the effective potential for the exciton 
relative motion quantum states with m r ^ has confinement energies larger than the one for the relative levels with 
m r = (compare, e.g., the energy range in Fig. 1(a) and Fig. 1(b)). 

The last term in our effective potential, Eq. (|20|), is related to the angular momentum of the exciton center-of-mass 
motion, tiir. In that case, it contributes to the potential like a centrifugal term (mu/R 2 ), and leads to a peak near 
R = (see dotted curve in Fig. 1-2 (b)). Also notice that the effective mass is position dependent (see dashed curve 
in Fig. l-2(b)). It is larger than the normal exciton mass M = 0.407 although its increase is small. 



B. Superconducting Disk 



The magnetic antidot was recently discussed by Reijniers et a£3, where the properties of electrons in a two di- 
mensional system confined by such a magnetic antidot were investigated. They found that, the magnetic field profile 
below a very thin type-I superconducting disk is given by: 



B z (R)=B a {l + - 
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with 
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where a is the disk radius, and d the distance of the superconducting disk to the 2DEG. In Fig. 3(a) the magnetic 
field profile resulting from the superconducting disk (magnetic antidot), and the effective potential, Eq. (pOJ), for the 
Is exciton relative motion quantum state is shown for a = 0.3 /jm and for d/a = 0.1 (dashed and dashed-dot curves) 
and d/a = 0.01 (solid and dotted curves). The effective mass, Eq. ([19]), and the effective potential for the 2p~ exciton 
relative motion quantum state is plotted in Fig. 3(b) for d/a = 0.1. We took in Fig. 3(a) (Fig. 3(b)) an applied 
magnetic field B a of 2 T (1 T). Notice that the magnetic field in the 2DEG is suppressed below the superconducting 
disk and an overshoot is found near the edge of the disk which becomes smoother with increasing distance, d, between 
the superconducting disk and the 2DEG. The exciton will prefer to localize below the disk because V e ^ is smaller 
there. This is no longer true for the 2p~ state (see Fig. 3(b)) where now V e 'f has a local minimum near the edge 
of the disk where the magnetic field exhibits the largest gradient. This occurs when the applied field B a is such that 
the diamagnetic term (quadratic in the magnetic field) in the effective potential, Eq. (|20|), becomes comparable to 
the angular momentum term (linear in the magnetic field). 



V. NUMERICAL RESULTS FOR THE EXCITON TRAPPING ENERGY 



We have calculated the trapping energy and wave functions of an exciton in a GaAs 2DEG in the presence of the 
above non-homogeneous magnetic field profiles (magnetic dipole and magnetic antidot profiles). The electron and 
hole mass and the dieletric constant used in our calculations area m e = 0.067mo, rrih = 0.34mo, and e — 12.53, 
respectively. The problem of finding the eigenfunctions and energies of Eq. ( |l8| ) was solved by using a similar 
numerical discretization technique as was used in Ref. 20. 

We define the exciton trapping energy as the difference between the exciton energy in the homogeneous applied field 
B a and the energy of the exciton in the non-homogeneous magnetic field profile for the same exciton state. Following 
this definition, a positive exciton trapping energy implies that the exciton is trapped in the effective potential created 
by the field inhomogeneity. The influence of all the terms in Eqs. (21) and (22) on the exciton trapping are extensively 
discussed. The possibility of exciton trapping by using the confinement potential created by the magnetic dipole, Figs. 
1-2, and the magnetic antidot, Fig. 3, are analyzed. In all the following figures, if not explicitly stated otherwise, 
the disk radius and the distance to the 2DEG used for the magnetized disk and for the superconducting disk are 
a = 2 fim, and a = 0.3 /xm, respectively, and we took the distance to the 2DEG equal to d = 0.1a. 

The trapping energy for the ground state of the exciton, for the center-of-mass quantum numbers (hr, tur) = (1,0) 
(dashed curves) and (1,1) (triangles) are shown in Fig. 4 for the magnetic dipole profile with a homogeneous applied 
field of (a) B a = 0.35 T and (b) B a = — 0.25T as a function of the magnetisation of the disk Bq . Similar results 
are shown in Fig. 5 for the magnetic antidot profile as a function of the applied magnetic field B a . Notice that 
for both profiles in low fields, the trapping energy is negative (i.e., unbound state), indicating that a large part of 
the center-of-mass exciton wave function is extended into the magnetic barrier region. The (nn,mji) = (1,1) is an 
excited state and requires slightly larger magnetic fields to become bound. With increasing magnetic field the exciton 
becomes more and more confined which increases the trapping energy Eb- The magnetic disk in the presence of a 
negative applied field (Fig. 4(b)) has a critical field where the trapping energy starts to decrease with increasing 
field. This is due to the competition between the magnetic field generated by the disk Bq , and the applied field B a . 
When the field of the disk increases, the center region of the corresponding effective potential (see solid curve in Fig. 
2(a)) which has a peak structure can become comparable to y e " at the edge of the disk which for large Bq (e.g. 
for (njj, ran) = (1,0) in Fig. 4(b) this occurs for Bq > 0.06 T) may lead to a complete vanishing of the confinement 



region. One can increase the trapping energy by applying a stronger homogeneous field B ai which is not trivial due 
to the limited coercivity field and for sufficient large B a it will flip the magnetisation of the disk and we end up into 
the situation of Fig. 1(a). 

The angular momentum interaction with .-the magnetic field is responsible for important effects in the trapping of 
atoms by non-homogeneous magnetic fieldsQ. In two-dimensional exciton systems, the momentum-field interaction 
is of importance in a variety of exciton properties. The exciton trapping energies for the seven lowest levels of the 
center-of-mass motion, for the (a) 2p~, (b) 2s, and (c) 2p + , exciton relative motion quantum states are shown in 
Fig. 6 (Fig. 7) for the magnetic dipole profile with B a = 0.35 T (B a = — 0.25T) as a function of Bq , and in Fig. 
8 for the magnetic antidot as a function of B a . In the case of the magnetic dipole profile, the exciton in the 2p ± 



states are much more confined than the 2s state. The reason is that the angular momentum term in Eq. (|20|) gives a 
confinement potential stronger than the diamagnetic term. Notice that in the case of a negative applied field, due to 
this strong confinement of the 2p ± states (Fig. 7), only the 2s state is affected by the competition between the fields 
Bq and B a as was discussed for Fig. 3. The trapping energy of the superconducting disk profile also increases with 
increasing magnetic field but for the 2p~ exciton relative motion quantum state (Fig. 8(a)) there is a local maximum 
in the trapping energy. This energy starts to decrease after some field (B a = 1 T in Fig. 8(a)), which is exactly when 
the quadratic term (positive) begins to be comparable to the angular momentum term (negative) in Eq. (EM (see 
effective potential in Fig. 3(b)). 

Also notice that the exciton center-of-mass quantum levels with non zero angular momentum (m,R ^ 0) have only 
slightly higher energy than the levels with m^j = 0. The reason is because the rriR — term in the effective potential, 
Eq. (EQ), only is responsible for a very small peak near R = which gives a very small contribution to the total 
effective potential. The trapping energy in Figs. 6-8 are typically two orders of magnitude larger than the one for 
the Is state of the exciton relative motion. In general, the trapping energy corresponding to the n r exciton relative 
quantum state resulting from the magnetic field inhomogeneity should be bigger than the one corresponding to the 
n r — 1 quantum state. This can be easily understood from the effective potential equation, Eq. (pH). The coefficient 
of the diamagnetic term, /3^ , increases with increasing n r which increases the exciton trapping energy. 

Figs. 9(a,b) and Fig. 10 show the ground-state radial wave function of the exciton center-of-mass motion (see Eq. 
(O)) for the magnetic dipole profile for B a — 0.35 T and B a = —0.25 T, and for the superconducting disk, respectively, 
as a function of the radial coordinate R, for the quantum numbers of the relative motion (n r , m r ) = (1,0), for different 
values of B® (magnetized disk) and B a (superconducting disk). Notice that with increasing magnetic field the exciton 
becomes more localized. The wave functions in Figs. 9(b) and 10 correspond to the effective potential profiles of Figs. 
2(a) and Fig. 3(a) which have the minimum of the potential near the center of the disk, and that is where the exciton 
is localized. In contrast to the case of Fig. 9(a) the effective potential (see Fig. 1(a)) has its minimum near the outer 
edge of the disk where the exciton becomes localized. In the latter case the center-of-mass wave function will have a 
ring like structure. 

In order to investigate the dependence of the exciton trapping energy on the size of the magnetic field inhomogeneity, 
we calculated the trapping energy for the 2s exciton relative motion quantum state (i.e., (n r , m r ) = (2,0)), as a function 
of the disk radius a. The results are shown in Figs. ll(a,b) for the magnetic dipole profile with Bq = 0.05 T and 
a homogeneous field of B a — 0.35 T and B a — — 0.25T, respectively, and in Fig. 12 for the magnetic antidot profile 
with Bq = 0.5 T, in each case for the seven lowest quantum numbers of the center-of-mass motion. In Fig. 11 (Fig. 
12) we took the distance to the 2DEG d = 0.2 /im (0.03 fim). The trapping energy increases with disk radius but it 
saturates for large a. The reason is that for sufficient large a only the width of the potential minimum increases but 
no longer the depth and consequently the trapping energy becomes equal to the depth of the minimum in V e *'. 

VI. CONCLUSIONS 

In conclusion, we have developed a formalism to describe the exciton motion in a two dimensional system, in the 
presence of a non-homogeneous magnetic field. We assumed that the length scale for the variation of the magnetic field 
is large as compared to the exciton radius. We have theoretically calculated the exciton trapping energies for the non- 
homogeneous magnetic fields created by a magnetized disk (the magnetic dipole type profile) and by a superconducting 
disk (the magnetic antidot). The results shown than the trapping energy has a significant dependence on the radius 
and profile of the non-homogeneous magnetic field, as well as on the homogeneous applied background magnetic field 
B a . Further, the exciton angular momentum interaction with the non-homogeneous magnetic field can be responsible 
for shifts as big as a factor of 10 in comparison with the states with zero angular momentum. In the present analysis 
we neglected the spin of the electron and the hole which may have an important effect on the trapping energy when 
the effective (/-factor is substantially different from zero. The analysis of the latter will be left for future research. 
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FIG. 1. (a) Effective potential for mn — and (n r ,m r ) = (1,0) and the magnetic field profile B(R) of the magnetized 
disk (inset of figure magnetic field profile in the 2DEG in the absence of a background field) with a positive homogeneous 
applied field B a — 0.35 T, as a function of the radial coordinate R, for d/a = 0.1 with V"'* (R) (solid) and B(R) (dotted), and 
d/a = 0.2 with V eff (R) (dashed-dotted) and B(R) (dashed), (b) Effective potential V eff (R) for n r = 2, m r = -1 with: i) 
ran = (solid) with corresponding effective mass M e ^(R) (dashed), and ii) V e ^(R) for tub, — 1 (dotted) as a function of R, 
for d/a — 0.1. We took Bq — 0.1 T and the disk radius a — 2.0 /an. 

FIG. 2. (a) Effective potential for mu = and (n r ,m r ) = (1,0) and the magnetic field profile B(R) of the magnetized disk 
with a negative homogeneous applied field B a = — 0.25T, as a function of the radial coordinate R, for d/a — 0.1 with V e '*{R) 
(solid) and B(R) (dotted), and d/a = 0.2 with V eff (R) (dashed-dotted) and B(R) (dashed), (b) Effective potential V eff (R) 
for n r = 2, m r = — 1 with: i) ran — (solid) with corresponding effective mass M e '' (R) (dashed), and ii) V £ '' (R) for vhr = 1 
(dotted) as a function of R, for d/a = 0.1. We took B® = 0.1 T and the radius a = 2.0 fim. 



FIG. 3. (a) Effective potential V e ^(R) for tur = 0, (n r ,m r ) — (1,0) and corresponding magnetic field profile B(R) of the 
superconducting disk, as a function of the radial coordinate 7?, for d/a — 0.1 with V e ** (R) (dashed-dot) and B(R) (dashed), and 
d/a = 0.01 with V eff (R) (solid) and B(R) (dotted), and (b) effective potential V eff (R) for n r = 2, m r = -1, m R = (solid) 
with corresponding effective mass M e ^ (R) (dashed), and V e ^ (R) for uir = 1 (dotted), as a function of R, for d/a — 0.1. We 
took the radius a — 0.3 /im and an applied field Bo of (a) 2 T and (b) 1 T. 

FIG. 4. Exciton trapping energy, Eb, for (nn,mR) = (1,0) (dashed curve) and (1,1) (triangles), and (n r ,m r ) — (1,0), for 
the magnetic dipole profile as a function of Bq , for (a) B a — 0.35 T and (b) B a — — 0.25T, and for a = 2 pm and d/a — 0.1. 

FIG. 5. Exciton trapping energy, Eb, for (nR,mR) = (1,0) (dashed curve) and (1,1) (triangles), and (n r ,m r ) — (1,0), for 
the magnetic antidot profile as a function of B a , and for d/a = 0.1, with a — 0.3 /im. 

FIG. 6. Exciton trapping energy, Eb, as a function of B$ , for the magnetic dipole profile for a homogeneous applied field 
B a — 0.35 T, for a = 2/im and d/a = 0.1, and for n r — 2, and (a) m r = —1, (b) m r = 0, and (c) m r — +1. The different 
curves correspond to the different exciton center-of-mass levels {nR,mn) = (1,0) dashed curve, (1,1) triangles, (2,0) dotted 
curve, (2,1) squares, (3,0) solid curve, (3,1) circles, and (4,0) dashed-dotted curve. 

FIG. 7. The same as in Fig. 6 but now for the magnetic dipole profile with a homogeneous applied field of B a — —0.25 T. 

FIG. 8. The same as in Fig. 6 but now as a function of B a , for the magnetic antidot profile, for a = 0.3 /im and d/a — 0.1. 

FIG. 9. Exciton center-of-mass radial wave function ip(R) as a function of the radial coordinate R, for the magnetic dipole 
profile, for (a) B a = 0.35 T and (b) B a — — 0.25T, for (n,R,mii) = (1,0), (n r ,m r ) = (1,0), for a — 2/im and d/a — 0.1, and for 
B ( f = 0.02 T (dotted), 0.06 T (solid), and 0.1 T (dashed). 

FIG. 10. The same as in Fig. 9(a) but now for the magnetic antidot profile, and for a = 0.3 /im and d/a = 0.1, and for 
B a = lT (dotted), 1.5 T (solid), and 2.0 T (dashed). 

FIG. 11. Exciton trapping energy, Eb, for (n r ,m r ) = (2,0), as a function of the disk radius a, for the magnetic dipole profile 
for (a) B a = 0.35 T and (b) B a — —0.25 T, for d = 0.2 /im and By — 0.05 T, and for center-of-mass quantum numbers (n_R, iur) 
= (1,0) dashed curve, (1,1) triangles, (2,0) dotted curve, (2,1) squares, (3,0) solid curve, (3,1) circles, and (4,0) dashed-dotted 
curve. 



FIG. 12. The same as in Fig. 11(a) but now for the magnetic antidot profile, and for d = 0.03 /im and B a — 0.5 T. 
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